Topological phase separation in 2D quantum lattice Bose-Hubbard system away from 

half-filling 
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We suppose that the doping of the 2D hard-core boson system away from half-filling may result in 
the formation of multi-center topological inhomogeneity (defect) such as charge order (CO) bubble 
domain(s) with Bose superfluid (BS) and extra bosons both localized in domain wall(s), or a topo- 
logical CO+BS phase separation, rather than an uniform mixed CO+BS supersolid phase. Starting 
from the classical model we predict the properties of the respective quantum system. The long- 
wavelength behavior of the system is believed to remind that of granular superconductors, CDW 
materials, Wigner crystals, and multi-skyrmion system akin in a quantum Hall ferromagnetic state 
of a 2D electron gas. 

To elucidate the role played by quantum effects and that of the lattice discreteness we have 
addressed the simplest nanoscopic counterpart of the bubble domain in a checkerboard CO phase of 
2D hc-BH square lattice. It is shown that the relative magnitude and symmetry of multi-component 
order parameter are mainly determined by the sign of the nn and nnn transfer integrals. In general, 
the topologically inhomogeneous phase of the hc-BH system away from the half-filling can exhibit 
the signatures both of s, d, and p symmetry of the off-diagonal order. 



I. INTRODUCTION 



The model of quantum lattice Bose gas has a long his- 
tory and has been suggested initially for conventional 
superconductors^, quantum crystals such as 4 He where 
superfluidity coexists with a crystalline orderiS^ After- 
wards, the Bose-Hubbard (BH) model has been stud- 
ied as a model of the superconductor-insulator transi- 
tion in materials with local bosons, bipolarons, or pre- 
formed Cooper pairsA^ Two-dimensional BH models 
have been addressed as relevant to describe the super- 
conducting films and Josephson junction arrays. The 
most recent interest to the system of hard-core bosons 
comes from the delightful results on Bosc-Einstcin (BE) 
condensed atomic systems produced by trapping bosonic 
neutral atoms in an optical latticed The progress in bo- 
son physics generates an especial interest around nonlin- 
ear topological excitations (skyrmions, vortices) which 
play an increasingly significant role in physics. 

One of the fundamental hot debated problems in 
bosonic physics concerns the evolution of the charge or- 
dered (CO) ground state of 2D hard-core BH model 
(hc-BH) with a doping away from half-filling. Numer- 
ous model studies steadily confirmed the emergence of 
"supersolid" phases with simultaneous diagonal and off- 
diagonal long range order while Penrose and Onsager— 
were the first showing as early as 1956 that supersolid 
phases do not occur. The most recent quantum Monte- 
Carlo (QMC) simulations&SiiS found two significant fea- 
tures of the 2D Bose-Hubbard model with a screened 
Coulomb repulsion: the absence of supersolid phase at 
half-filling, and a growing tendency to phase separation 
(CO+BS) upon doping away from half-filling. Moreover, 
Batrouni and Scalettas^ studied quantum phase transi- 
tions in the ground state of the 2D hard-core boson Hub- 
bard Hamiltonian and have shown numerically that, con- 



trary to the generally held belief, the most commonly dis- 
cussed "checkerboard" supersolid is thermodynamically 
unstable and phase separates into solid and superfluid 
phases. The physics of the CO+BS phase separation in 
Bose-Hubbard model is associated with a rapid increase 
of the energy of a homogeneous CO state with doping 
away from half-filling due to a large " pseudo-spin-flip" 
energy cost. Hence, it appears to be energetically more 
favorable to " extract" extra bosons (holes) from the CO 
state and arrange them into finite clusters with a rela- 
tively small number of particles. Such a droplet scenario 
is believed to minimize the long-range Coulomb repul- 
sion. However, immediately there arise several questions. 
Whether a simple mean-field approximation (MFA) and 
classical continuum model can predict such a behavior? 
What is the detailed structure of the CO+BS phase sep- 
arated state? What are the main factors governing the 
essential low-energy and long-wavelength physics? Is it 
possible to make use of simple toy models? Are there 
the analogies with fermion Hubbard model? The behav- 
ior of the latter under doping away from half-filling is 
extensively studied in last decade in frames of stripe sce- 
nario which implies that doping may proceed through 
the formation of stripes, or charged domain walls being 
specific topological solitons in Neel antiferromagnets. In 
particular, Emery and Kivelson 11 argued that the phase 
separation reflects a universal tendency of the correlated 
antifcrromagnet to expel the doped holes. White and 
Scalapinoi^ showed that the pure 2D t-J model, in the 
small- J/t regime and with dopings near x ~ 0.1, has a 
striped ground state. In general, the topological phase 
separation is believed to be a generic property of 2D 
fermion Hubbard model. 

In the paper we present a topological scenario of 
CO+BS phase separation in 2D hc-BH model with inter- 
site repulsion. The extra bosons/holes doped to a 
checkerboard CO phase of 2D boson system are believed 
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to be confined in the ring-shaped domain wall of the 
skyrmion-likc topological defect which looks like a bub- 
ble domain in an easy- axis antifcrromagnet. This allows 
us to propose the mechanism of 2D topological CO+BS 
phase separation when the doping of the bare checker- 
board CO phase results in the formation of a multi-center 
topological defect, which simplified pseudo-spin pattern 
looks like a system of bubble CO domains with Bose su- 
pcrfluid confined in charged ring-shaped domain walls. 

The rest of the paper is organized as follows. In Sec. II 
we give a short overview of the conventional hc-BH model 
in frames of a pseudo-spin formalism and MFA descrip- 
tion. In Sec. Ill we show that the doping in 2D hc-BH 
system can be accompanied by the formation of a topo- 
logical defect like a single or multi-center skyrmion. We 
present the scenario of the essential low-energy physics 
for the topologically doped hc-BH system. In Sec. IV we 
address a simple model of nanoscopic bubble-like domain 
in a checkerboard CO phase for a discrete square lat- 
tice that allows us to demonstrate a subtle microscopic 
structure of such a center with a multi-component order 
parameter. 

II. HARD-CORE BOSE-HUBBARD MODEL 
A. Effective Hamiltonian. Pseudo-spin formalism 

The Hamiltonian of hard-core Bose gas on a lattice can 
be written in a standard form as follows: 

h bg Y^i.J'uku, ■ //;/*, = /' -^]r,, .v. .v, i>y,\,. 

i>j i>j i 

(1) 

where P is the projection operator which removes dou- 
ble occupancy of any site, JVf = Sj-Bi, /i the chem- 
ical potential determined from the condition of fixed 

N 

full number of bosons Ni = {Ni) or concentration 

i=l 

n = Ni/N e [0,1]. The tij denotes an effective trans- 
fer integral, is an intersite interaction between the 
bosons. Here B^(B) are the Pauli creation (annihila- 
tion) operators which are Bose-like commuting for dif- 
ferent sites [Bi,Bt] = 0, if i 7^ j, [B^Bj] = 1 - 2N l: 

Ni = BjBi; N is a full number of sites. It is worth not- 
ing that near half- filling (n k, 1/2) one might introduce 
the renormalization Ni — > (Ni — 1/2), or neutralizing 
background, that immediately provides the particle-hole 
symmetry. 

The model of hard-core Bose-gas with an intersite re- 
pulsion is equivalent to a system of spins s = 1/2 exposed 
to an external magnetic field in the z-direction. For the 
system with neutralizing background we arrive at an ef- 
fective pseudo-spin Hamiltonian 

Hbg = J??(4sj + afar) + g j* - ^ ** , 

i>j i>j i 

(2) 



where J%> = 2ty, J?. = Vy, s~ = ±Bs+ = 

-rk B ^ sZ = -i + B \ B ^ s± = T71 (* x ± **")■ 

In a linear approximation the Hamiltonian for the cou- 
pling with an electromagnetic field reads as follows 

vm Xa.- ::<|, < - **)[^ x S ^ - - $ *) 2 ( § * • s j)); 

i<j 

(3) 

($ j -$ i ) = --jLj* 3 A(r)dr, (4) 

where A is the vector-potential, and integration runs over 
line binding the i and j sites; 

(Si -Si) = \{ B \ B 3+ B i B ]), [sixsjl = l -(BjB 3 -B t Bf), 

(5) 

however, the pseudo-spins are assumed to lie in xy plane. 
Then the boson current density operator one may repre- 
sent to be a sum of the paramagnetic 

y(Ri) = I ^t«Rtf[si x s j]*> ( 6 ) 

3 

and diamagnetic 

terms, respectively. 

B. Uniform phases: mean-field approximation 

Below we make use of a conventional two-sublattice 
MFA approach. For the description of the pseudospin or- 
dering to be more physically clear one may introduce two 
vectors, the ferromagnetic and antiferromagnetic ones: 

m = ^«si) + (sa)); 1 = -^((si) - (s 2 »; m 2 + l 2 = 1 . 

where m • 1 = 0. For the plane where these vec- 
tors lie one may introduce two-parametric angular de- 
scription: m x = sin a cos/3, m z = — sin a sin (3, l x = 
cos a sin /3, l z = cos a cos (3. The hard-core boson system 
in a two-sublattice approximation is described by two 
diagonal order parameters l zi m z and two complex off- 
diagonal order parameters m± — ^f^(m x ± im y ) and 

l± = ^f-±=(l x ± il y ). The complex superfluid order pa- 
rameter ^(r) = |^(r)| exp(— up) is determined by the 
in-plane components of ferromagnetic vector: ^(r) = 
\((B\ + B 2 )) = sni- = sm±_ exp(— up), m± being the 
length of the in-plane component of ferromagnetic vec- 
tor. So, for a local condensate density we get n s = s 2 m 2 L . 
It is of interest to note that in fact all the conventional 
uniform T = states with nonzero ^(r) imply simul- 
taneous long-range order both for modulus |^(r)| and 
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phase ip. The in-plane components of antiferromagnetic 
vector l± describe a staggered off-diagonal order. It is 
worth noting that by default one usually considers the 
negative sign of the transfer integral , that implies the 
ferromagnetic in-plane pseudospin ordering. 

The model exhibits many fascinating quantum phases 
and phase transitions. Early investigations predict at 
T = charge order (CO), Bose superfluid (BS) and 
mixed (CO+BS) supersolid uniform phases with an 
Ising-type melting transition (CO-NO) and Kosterlitz- 
Thouless-type (BS-NO) phase transitions to a non- 
ordered normal fluid (NO)^& The detailed mean-field and 
spin-wave analysis of the uniform phases of 2D hc-BH 
model is given by Pich and Freyi 1 ^ MFA yields for the 
conventional uniform supersolid phase 4 



• 2 fl VV^Zi . 2 
sin p = m z t= =; sin a 



WV + 2t : 



m. 



yv + 2t 

WV-2t 



with a constant chemical potential fi = AsyJ{V 2 — At 2 ). 
It should be noted that the supersolid phase appears to 
be unstable with regard to small perturbations in the 
Hamiltonian. The mean-field energy per site of the uni- 
form supersolid phase is written as follows: 

Ess = E C o + s[im z = E C o + K n B ~ \), 

where E C o = —2Vs 2 . The cost of doping both for CO 
and CO+BS phase appears to be rather high as com- 
pared with an easy-plane BS phase at half-filling where 
the chemical potential turns into zeroi^ 



III. DOPING OF BH SYSTEM AWAY FROM 
HALF-FILLING: CONTINUOUS LIMIT 



Topological phase separation: skyrmion-like 
bubble domains 



Magnetic analogy allows us to make unambiguous pre- 
dictions as regards the doping of BH system away from 
half-filling. Indeed, the boson/hole doping of checker- 
board CO phase corresponds to the magnetization of an 
antifcrromagnet in z-direction. In the uniform easy-axis 
Z 2 -phase of anisotropic antiferromagnet the local spin- 
flip energy cost is rather big. In other words, the en- 
ergy cost for boson/hole doping into checkerboard CO 
phase appears to be big due to a large contribution of 
boson-boson repulsion. However, the magnetization of 
the anisotropic antiferromagnet in an easy axis direc- 
tion may proceed as a first order phase transition with a 
"topological phase separation" due to the existence of an- 
tiphase domains. 16 The antiphase domain walls provide 
the natural nucleation centers for a spin-flop phase hav- 
ing enhanced magnetic susceptibility as compared with 
small if any longitudinal susceptibility thus providing 
the advantage of the field energy. Namely domain walls 
would specify the inhomogeneous magnetization pattern 



for such an anisotropic easy-axis antiferromagnet in rel- 
atively weak external magnetic field. As concerns the 
domain type in quasi-two-dimensional antiferromagnet 
one should emphasize the specific role played by the 
cylindrical, or bubble domains which have finite energy 
and size. These topological solitons have the vortex-like 
in-plane spin structure and often named " skyrmions" . 
The classical, or Belavin-Polyakov (BP) skyrmionii de- 
scribes the solutions of a non-linear cr-model with a clas- 
sical isotropic 2D Heisenberg Hamiltonian and represents 
one of the generic toy model spin textures. It is of pri- 
mary importance to note that skyrmion-like bubble do- 
mains in easy-axis layered antiferromagnets were actually 
observed in the experiments performed by Waldner, 18 
that were supported later by different authors (see, e.g. 
Refs ll9l20l) . Although some questions were not com- 
pletely clarified and remain open until nowf2L2£ the clas- 
sical and quantum2ii2i skyrmion-like topological defects 
are believed to be a genuine element of essential physics 
both of ferro- and antiferromagnetic 2D easy-axis sys- 
tems. 

The magnetic analogy seems to be a little bit naive, 
however, it catches the essential physics of doping the 
hc-BH system. As regards the checkerboard CO phase of 
such a system, namely a finite energy skyrmion-like bub- 
ble domain^ 5 - seems to be the most preferable candidate 
for the domain with antiphase domain wall providing the 
natural reservoir for extra bosons. The classical descrip- 
tion of nonlinear excitations in hc-BH model includes the 
skyrmionic solution given V = 2tm^ The skyrmion spin 
texture consists of a vortex-like arrangement of the in- 
plane components of ferromagnetic m vector with the l z - 
component reversed in the centre of the skyrmion and 
gradually increasing to match the homogeneous back- 
ground at infinity. The simplest spin distribution within 
classical skyrmion is given as follows 

m x — m± cos(f + ipo); m y = m± sin(y + <po); 



2r\ 
r 2 + X 2 1 



h = 



A 2 



(8) 



where cpo is a global phase (t/(l) order parameter), A 
skyrmion radius. The skyrmion spin texture describes 
the coexistence and competition of the staggered charge 
order parameter l z and BS order parameter m± (see 
FigQ that reflects a complex spatial interplay of poten- 
tial and kinetic energies. The skyrmion looks like a bub- 
ble domain in an easy-axis magnet. It should be noted 
that the domain wall in such a bubble domain somehow 
created in the checkerboard CO phase of 2D hc-BH sys- 
tem represents an effective ring-shaped reservoir both for 
Bose-superfluid and extra boson/hole. Indeed, the soli- 
ton energy depends quadratically on the number An of 
bosons bound in domain walljSSi similarly to that of ho- 
mogeneous BS phase. 15 In other words, one might say 
about a zero value of the effective boson/hole chemical 
potential for the CO bubble domain configuration, pro- 
vided it were a ground state. The numerical calculations 
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FIG. 1: The order parameters distribution in the skyrmion: 
a) The radial distribution of the staggered charge order pa- 
rameter l z ; b) the radial distribution of the modulus of the 
superfluid order parameter mx; c) the radial distribution of 
the charge density m z for the charged skyrmion; d) the orien- 
tation of the superfluid order parameter mx given two values 
of the phase order parameter tpo: tpo — 0,7r/2, respectively. 
Rings in all the pictures correspond to the skyrmion (bubble 
domain) wall, the dashed circle being its center. 



performed in frames of a classical continuous models 
show that the doped bosons appear to be trapped in- 
side the bubble domain wall. The spatial distribution 
of the doped boson/hole density (oc m z ) in the charged 
skyrmion is shown schematically in Fig. lc. 

Skyrmions are characterized by the magnitude and 
sign of its topological charge, by its size (radius), and 
by the global orientation of the spin. The scale invari- 
ance of classical BP skyrmionic solution reflects in that 
its energy does not depend on radius, and global phase. 
An interesting example of topological inhomogeneity is 
provided by a multi-center BP skyrmionii which energy 
does not depend on the position of the centers. The lat- 
ter are believed to be addressed as an additional degree 
of freedom, or positional order parameter. 

In the continuous model the classical BP skyrmion is a 
topological excitation and cannot dissipate. However, the 
classical static skyrmion is unstable with regard to an ex- 
ternal field or anisotropic interactions both of easy-plane 
and easy-axis type. Small easy-axis anisotropy or exter- 
nal field are sufficient to shrink skyrmion to a nanoscopic 



size when magnetic length 

z = (vW*) 2 - w*) 2 - 4 ) 

is of the order of several lattice parameters, and the 
continuous approximation fails to correctly describe ex- 
citations. Nonetheless, Abanov and PokrovskySi have 
shown that the easy-axis anisotropy together with fourth- 
order exchange term can stabilize skyrmion with radius 
R oc 

B. Topological phase separation: Skyrmion lattices 
and the low-energy physics of BH model away from 
half-filling 

A skyrmionic scenario in hc-BH model allows us to 
make several important predictions. Away from half- 
filling one may anticipate the nucleation of a topologi- 
cal defect in the unconventional form of the multi-center 
skyrmion-like object with ring-shaped Bose superfluid re- 
gions positioned in an antiphase domain wall separating 
the CO core and CO outside of the single skyrmion. The 
specific spatial separation of BS and CO order parame- 
ters that avoid each other reflects the competition of ki- 
netic and potential energy. Such a topological (CO+BS) 
phase separation is believed to provide a minimization 
of the total energy as compared with its uniform super- 
solid counterpart. Thus, the parent checkerboard CO 
phase may gradually lose its stability under boson/hole 
doping, while a novel topological self-organized texture 
is believed to become stable. The most probable possi- 
bility is that every domain wall accumulates single bo- 
son, or boson hole. Then, the number of centers in a 
multi-center skyrmion nucleated with doping has to be 
equal to the number of bosons/holes. In such a case, 
we anticipate the near-linear dependence of the total BS 
volume fraction on the doping. Generally speaking, one 
may assume scenario when the nucleation of a multi- 
center skyrmion occurs spontaneously with no doping. In 
such a case we should anticipate the existence of neutral 
multi-center skyrmion-like object with equal number of 
positively and negatively charged single skyrmions. How- 
ever, in practice, namely the boson/hole doping is likely 
to be a physically clear driving force for a nucleation of a 
single, or multi-center skyrmion-like self-organized collec- 
tive mode in the form of multi-center charged topological 
defect which may be (not strictly correctly) referred to 
as multi-skyrmion system akin in a quantum Hall fer- 
romagnetic state of a two-dimensional electron gas^ In 
such a case, we may characterize an individual skyrmion 
by its position (i.e., the center of skyrmionic texture), 
its size (i.e., the radius of domain wall), and the orien- 
tation of the in-plane components of pseudo-spin (U(l) 
degree of freedom) . An isolated skyrmion is described by 
the inhomogeneous distribution of the CO parameter, or 
staggered boson density l z , order parameter m z charac- 
terizing the deviation from the half-filling, and m± that 



5 



corresponds only to the modulus of the superfluid order 
parameter. 

It seems likely that for a light doping any doped parti- 
cle (boson/holes) results in a nucleation of a new single- 
skyrmion state, hence its density changes gradually with 
particle doping. Therefore, as long as the separation be- 
tween skyrmionic centers is sufficiently large so that the 
inter-skyrmion interaction is negligible, the energy of the 
system per particle remains almost constant. This means 
that the chemical potential of a boson or hole remains un- 
changed with doping and hence apparently remains fixed. 

The multi-skyrmionic system in contrast with uniform 
ones can form the structures with inhomogeneous long- 
range ordering of the modulus of the superfluid order 
parameter accompanied by the non-ordered global phases 
of single skyrmions. Such a situation resembles in part 
that of granular superconductivity. 

In the long-wavelength limit the off-diagonal ordering 
can be described by an effective Hamiltonian in terms 
of U(l) (phase) degree of freedom associated with each 
skyrmion. Such a Hamiltonian contains a repulsive, long- 
range Coulomb part and a short-range contribution re- 
lated to the phase degree of freedom. The latter term can 
be written out in the standard for the XY model form 
of a so-called Josephson coupling 



Hj = - Jij cos(ip l -ipj), 



(9) 



where tpi, tpj are global phases for skyrmions centered at 
points i,j, respectively, Jij Josephson coupling param- 
eter. Namely the Josephson coupling gives rise to the 
long-range ordering of the phase of the superfluid order 
parameter in a multi-center skyrmion. Such a Hamil- 
tonian represents a starting point for the analysis of 
disordered superconductors, granular superconductivity, 
insulator-superconductor transition with array of 

superconducting islands with phases ipi,(pj. Calculating 
the phase-dependent part of skyrmion-skyrmion interac- 
tion Timm et alm^ arrived at negative sign of Jij that 
favors antiparallel alignment of the U(l) pseudospins. In 
other words, two skyrmions are believed to form a pe- 
culiar Josephson n micro-junction. There are a number 
of interesting implications that follow directly from this 
result:* 5 — the spontaneous breaking of time-reversal sym- 
metry with non-zero supercurrents and magnetic fluxes 
in the ground state, long-time tails in the dynamics of the 
system, unconventional Aharonov-Bohm period hc/Ae, 
negative magnetoresistence. 

To account for Coulomb interaction and allow for 
quantum corrections we should introduce into effective 
Hamiltonian the charging energy«2S 




where n% is a boson number operator for bosons bound 
in i-th skyrmion; it is quantum-mechanically conjugated 



bubble liquid with non- bubble crystal with 
ordered global phase global phase ordering 



FIG. 2: Bubble textures for the bubble liquid and bubble 
crystal. The arrows show the global phase order parameter. 
The left hand side panel shows snapshot of bubble texture in 
liquid state with non-ordered global phase. The right hand 
side panel illustrates the triangular bubble crystal state with 
" ferromagnetically" ordered global phase. Rings in all the 
pictures correspond to the bubble domain wall. 



to ip: m — —id/difi, (C 1 )y stands for the capacitance 
matrix, q for bosonic charge. 

The classical ground state energy of Skyrmion lat- 
tice for all reasonable two-dimensional lattice structures 
was minimized by Timm et al. taking the U(l) order 
into account^ Besides the expected triangular lattice 
with frustrated antiferromagnetic U(l) order and square 
lattice with Neel U(l) order the authors have also ob- 
tained ground-state energies for all 2D Bravais lattices. 
Such a system appears to reveal a tremendously rich 
quantum-critical structure. In the absence of disorder, 
the T — phase diagram of the multi-skyrmion system 
implies either triangular, or square crystalline arrange- 
ments (Skyrmion, or bubble crystal) with possible melt- 
ing transition to a Skyrmion (bubble) liquid. The melting 
of Skyrmion lattice is successfully described applying the 
Berezinsky-Kosterlitz-Thouless (BKT) theory to disloca- 
tions and disclinations of the lattice and proceeds in two 
steps. The first implies the transition to a liquid-crystal 
phase with short-range translational order, the second 
does the transition to isotropic liquid. Disorder pins the 
Skyrmion lattice and also causes the crystalline order to 
have a finite correlation length. For such a system pro- 
vided the skyrmion positions fixed at all temperatures, 
the long-wave-length physics would be described by an 
XY model with expectable BKT transition and gapless 
XY spin-wave mode. 

The classical phase diagram of Skyrmion lattice is quite 
rich. Depending on the relation between Coulomb and 
Josephson coupling, and density one may arrive at dif- 
ferent lattice structures with continuous or first order 
phase transitions. As regards the superfluid properties 
the skyrmionic system reveals unconventional behavior 
with two critical temperatures Tbs < t an d T c < J, Tbs 
being the temperature of the ordering of the modulus and 
T c < Tbs that of the phase of order parameter \P. The 
low temperature physics in Skyrmion crystals is governed 
by an interplay of two BKT transitions, for the U(l) 
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phase and positional degrees of freedom, respectively^ 
Dislocations in most Skyrmion lattice types lead to a mis- 
match in the U(l) degree of freedom, which makes the 
dislocations bind fractional vortices and lead to a cou- 
pling of translational and phase excitations. Both BKT 
temperatures either coincide (square lattice) or the melt- 
ing one is higher (triangular lattice)^ Quantum fluctu- 
ations can substantially affect these results. Quantum 
melting can destroy U(l) order at sufficiently low den- 
sities where the Josephson coupling becomes exponen- 
tially small. Similar situation is expected to take place 
in the vicinity of structural transitions in Skyrmion crys- 
tal. With increasing the skyrmion density the quantum 
effects result in a significant lowering of the melting tem- 
perature as compared with classical square-root depen- 
dence. The resulting melting temperature can reveal an 
oscilating behavior as a function of particle density with 
zeros at the critical (magic) densities associated with 
structural phase transitions. 

In terms of our model, the positional order corresponds 
to an incommensurate charge density wave, while the 
U(l) order does to a superconductivity. In other words, 
we arrive at a subtle interplay between two orders. The 
superconducting state evolves from a charge order with 
Tq < T m , where T m is the temperature of a melting 
transition which could be termed as a temperature of 
the opening of the insulating gap (pseudo-gap!?). 

The normal modes of a dilute skyrmion system (multi- 
center skyrmion) include the pseudo-spin waves propa- 
gating in-between the skyrmions, the positional fluctua- 
tions, or phonon modes of the skyrmions which are gap- 
less in a pure system, but gapped when the lattice is 
pinned, and, finally, fluctuations in the skyrmionic in- 
plane orientation and size. The latter two types of fluc- 
tuation are intimitely connected, since the z-component 
of spin and orientation are conjugate coordinates because 
of commutation relations of quantum angular momentum 
operators. So, rotating a skyrmion changes its size. The 
orientational fluctuations of the multi-skyrmion system 
are governed by the gapless XY model.— The relevant 
model description is most familiar as an effective theory 
of the Josephson junction array. An important feature of 
the model is that it displays a quantum-critical point. 

The low-energy collective excitations of skyrmion liq- 
uid includes an usual longitudinal acoustic phonon 
branch. The liquid crystal phases differ from the isotropic 
liquid in that they have massive topological excitations, 
i.e., the disclinations. One should note that the liquids 
do not support transverse modes, these could survive in 
a liquid state only as overdamped modes. So that it is 
reasonable to assume that solidification of the skyrmion 
lattice would be accompanied by a stabilization of trans- 
verse modes with its sharpening below melting transition. 
In other words an instability of transverse phonon modes 
signals the onset of melting. 

A generic property of the positionally ordered 
skyrmion configuration is the sliding mode which is usu- 
ally pinned by the disorder. The depinning of sliding 



mode(s) can be detected in a low-frequency and low- 
temperature optical response. 



C. Implications for the "doping-temperature" 
phase diagram of hc-BH model 

Our speculations as regards the topological phase sep- 
aration and long wave-length physics of 2D hc-BH model 
are summarized in a qualitative An-T (An;, = rif, — 0.5) 
phase diagram in Fig0 First of all, the phase diagram 
implies a scenario in which the topological phase sepa- 
ration state evolves with a minimal doping, though it 
is worth noting that the ideal checkerboard CO phase 
cannot be a nominal An;, = limit of any topologically 
phase separated phase. 

Despite the qualitative character of phase diagram we 
took into account some quantitative results of quantum 
Monte-Carlo calculations for the 2D hc-BH model with 
nearest neighbor Coulomb repulsion First, it con- 
cerns the doping range of additive, or near-linear concen- 
tration behavior of Bose-condensate density p s and CO 
structure factor 5(q). The superfluid density increases 
approximately linearly with doping except for the most 
overdoped point An& w 0.1, where it turns down again. 
Diagonal long-range charge order is characterized by the 
equal time structure factor at the ordering vector q, 

5(q) = ^E eiq Xj' T Hj + 1 ' T ))- ( 10 ) 
i 

In the presence of long-range order S*(q) will diverge with 
the system size <S(q) cx L 2 for a given ordering momen- 
tum Q, which characterizes the ordered phase. For the 
checkerboard order Q = (tt, 7t). In the concentration 
range Ant = 0.0 — 0.1, where both quantities vary lin- 
early with Aribfi we may approximate the topological 
defect to be a system of AA = NAni. interacting single- 
charged skyrmions. It should be noted that the both 
CO and BS order parameters coexist in a rather narrow 
doping concentration interval: An^ < 0.11. Beyond this 
" overdoping" region we deal with an inhomogeneous bo- 
son liquid which pseudo-spin picture implies a strongly 
frustrated singlet-triplet system that resembles the spin 
glass, and maybe termed as a dynamical "singlet-triplet" 
pseudo-spin glass. Such a system can be characterized 
by a dynamical short-range diagonal and off-diagonal or- 
dering with a wide distribution of respective correlation 
lengths and relaxation rates. Interestingly, that in frames 
of a continuous model this phase is still described to be a 
strongly correlated multi-center topological defect. How- 
ever, such a model fails obviously to describe the real 
system where the inter-center spacing is comparable with 
the lattice parameter. 

The temperature evolution of hc-BH system with large 
inter-site boson repulsion implies the highest critical tem- 
perature Tqo (An) separating the high-temperature non- 
ordered NO phase (boson liquid) and a low-temperature 
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FIG. 3: The low-temperature part of the schematical An-T phase diagram for the hc-BH model with topological phase 
separation. The gradual skyrmionic solidification evolves from the isotropic skyrmion(bubble) liquid phase, the skyrmion 
liquid-crystal phase, the incommensurate skyrmion crystal phase, and U(f) ordering, or superfluid phase, restricted by the 
temperatures of the proper BKT transitions. The right hand side panel differs by the assumed structural phase transition in 
the bubble lattice at a magic doping such as Ant = 1/16. 



quasi-CO phase with a disordered system of skyrmions. 
The next critical temperature Ttps (An) < t points to 
the first order phase transition with a formation of in- 
homogeneous Bose condensate in a single skyrmion with 
the vortex-like texture of the quasi-local order parameter 
mi. In other words, it is a temperature of the topolog- 
ical CO+BS phase separation (TPS). In frames of our 
scenario such TPS state emerges with a minimal doping, 
and TxpsiArib) is likely to be nearly constant in a linear 
doping regime (Ant, — 0.0 — 0.1). The low-temperature 
part of the phase diagram, which is schematically shown 
in Fig. |21 describes the gradual skyrmionic solidifica- 
tion and may include the isotropic skyrmion(bubble) liq- 
uid phase, the skyrmion liquid-crystal phase, the incom- 
mensurate skyrmion crystal phase, and U(l) ordering, 
or superfluid phase, restricted by the temperatures of 
the proper BKT transitions. For a small deviation from 
half- filling (" under-doping" ) the temperatures of bubble 
crystallization/melting are seemingly to obey the square 
root concentration law T rn oc V An& with a strongly de- 
veloped quantum melting effect when approaching an 
" overdoped" regime or concentrations limiting the linear 
regime. The superfluid phase in Figgis arbitrarily cho- 
sen to lie inside the skyrmion solid phase. One should 
emphasize the specific role played by quantum fluctua- 
tions: these lead to the melting of the bubble crystal at 
high densities and orientational disordering22i2i at low 
densities. Both effects are of primary importance in the 
overdoped and heavily underdoped regions of the phase 
diagram, respectively. Moreover, quantum melting effect 
may strongly affect the phase diagram near the "magic" 
doping level, where the skyrmion lattice undergoes the 



structural phase transition. For illustration, in the right 
hand side panel in Fig. |3we present the possible phase 
diagram for a hc-BH system with a quantum melting ef- 
fect near the "magic" doping level Ani, = 1/16. In all 
the cases, the critical doping level for the superfluid for- 
mation is determined by the magnitude of the Josephson 
coupling constant. It is worth noting that the bubble 
crystallization is accompanied by different (pseudo)gap 
effects. 

Our interpretation of the phase transition at Ttps dif- 
fers from that in Refllfl This temperature is governed by 
the magnitude of transfer integral, and believed to be a 
temperature of the emergence of the nonzero magnitude 
of the modulus of the superconducting order parameter 
rather than a critical temperature for an insulator to su- 
perconductor transition as it is stated in RefllfJ. This 
conclusion seems to be a result of finite size effects and 
boundary conditions in QMC calculations despite the au- 
thors make use of the most efficient QMC strategy. Such 
problems seem to be typical for a finite-size simulation of 
many phase transitions. In addition, we should note that 
MC calculations need substantial increase in lattice size 
to reproduce quantitatively the long-wavelength physics 
because the size of skyrmion and Skyrmion lattice pa- 
rameter appear to be new characteristic lengths. 

The close inspection of the phase diagram in FigEl 
where T c < T m < Ttps, does not provide the opti- 
mistic expectations as regards the high-temperature su- 
perconductivity in 2D hc-BH systems, even if the magni- 
tude of the local boson transfer integral were as large as 
t w 1000 K. Nevertheless, the attractively large temper- 
atures Ttps of the topological phase separation with the 
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emergence of nonzero local condensate density engender 
different reasonable speculations as regards its practical 
realization. 



IV. TOPOLOGICAL PHASE SEPARATION IN 
DISCRETE LATTICES 

The making use of the mean-field approximation to- 
gether with simplified classical continuous models can 
hardly provide the quantitative description of a quan- 
tum lattice boson system. Both quantum effects and the 
discreteness of skyrmion texture can result in substan- 
tial deviations from the predictions of classical model. 
The continuous model is relevant for discrete lattices 
only if we deal with long-wave length inhomogeneities 
when their size is much bigger than the lattice spac- 
ing. In the discrete lattice the very notion of topolog- 
ical excitation seems to be inconsistent. At the same 
time, the discreteness of the lattice itself does not pro- 
hibit from considering the nanoscale (pseudo)spin tex- 
tures whose topology and spin arrangement is that of a 
skyrmion^ Naturally, the typical continuous models fail 
to describe properly such short-wave length nanoscopic 
inhomogeneities. Hereafter, we discuss a simple model 
which seems to catch the main effects of discreteness and 
quantization. 



A. Nanoscopic bubble domain in checkerboard CO 

phase 

What is the lattice counterpart of the small skyrmion- 
like bubble domain? Fig. 01 presents the schematic 
view of such a smallest skyrmion-like bubble domain in 
a checkerboard CO phase for 2D square lattice with an 
effective size of 4 lattice spacings. The domain wall is 
believed to include as a minimum eight sites forming a 
ring-shaped system of four dimers each composed of two 
sites. There are two types of such domains which differ by 
a rotation by ±7r/2. The formation of such a center seems 
not to require a big energy. Indeed, the simple estimate 
of a change in potential energy: AV « 4 V nn — 6 V nnn 
points to a near cancelation of nn and nnn contribu- 
tions. The remarkable feature of the domain is in a rather 
small magnitude both of scalar potential and electric field 
inside the 8-site domain wall (see Fig. 0J. The flip- 
energy for the dimer dipole moment is estimated to be 
SV ~ V nn — 2 V nnn , compared with that for bare checker- 
board CO phase (SVq « 3V nn — 4 V nnn ) that implies 
rather subtle competition between nn and nnn couplings. 
However, the difference {SVq — SV) = 2{V nn — V nnn ) is be- 
lieved to be always positive and large. In other words, the 
flip-energy for the dimer dipole may be relatively small. 
Respective dipole fluctuations would result in an effective 
screening of electrostatic repulsion energy thus provid- 
ing a stabilization of a bubble- like defect. An important 
additional mechanism of the domain stabilization in the 




FIG. 4: Illustration to a skyrmion-like bubble domain in 
checkerboard CO phase with 8-site ring-shaped domain wall. 
The four dimers within the domain wall are marked by I-IV. 
Schematically shown are nearest and next-nearest neighbor 
domains which do not overlap with the central domain. It 
is worth noting that there are eight additional domains with 
strong nearest neighbor inter-dimer coupling with the central 
domain. The centers of all 12 nearest neighbor domains which 
do not overlap with the central domain are marked by crosses. 



hc-BH systems with local bosons composed from electron 
pairs may arise from the electron and lattice polarization 
effects M£l As is well known, the respective energies are 
comparable in value with the intersite Coulomb interac- 
tion. Probably, namely both these effects might strongly 
contribute to the domain stabilization energy. On the 
other hand, the bubble geometry implies the formation 
of the electrostatic potential well inside the domain wall 
both for positive and negative charges. It means that 
the doping into a domain wall stabilizes the domain con- 
figuration. Such a doping maybe energy costless, while 
the energy cost of the pseudo-spin-flip in a checkerboard 
CO phase is rather high: AV w 4V r „„ — 4V r „„„, if the 
neutralizing background is taken into account. 



B. Pseudo-spin formalism in a two-center dimer 

With taking into account the kinetic energy (quantum 
tunneling) only inside the 8-site domain wall we shall con- 
sider it as a quantum system in an external electrostatic 
potential field assuming a rigid checkerboard CO order- 
ing overall beyond the domain wall. Further, taking into 



account the reasonable relation: \t nn \ ^> \t r , 



we shall 



consider the domain wall to be a system of four dimers, 
or pairs of nearest neighboring sites forming a quantum 



9 



S=0 



S=l 



S=0 




isotropic! 
term 



lelectric ] 
Ifield (hj) 



effective 
magnetic 
field (h s ) 



FIG. 5: The step-by-step formation of the dimer energy spec- 
trum given different signs of the nn transfer integral. The 
arrows mark the dip ole- allowed transitions (see text for de- 
tails). 



cluster. We anticipate the model though being simplified 
would be very instructive analyzing the role played both 
by quantum and discreteness effects. 

Let first address a simple two-site, or a dimer system. 
The effective pseudo-spin-Hamiltonian for such a cluster, 
or a single dimer Hamiltonian can be written as follows 



H d = -t nn [S(S- 



DS 2 Z -h s S z - h T f z , (11) 



where S = §i + §2 is the total pseudo-spin momentum, 
T = §1 — §2 is an operator that changes pseudo-spin mul- 
tiplicity, D = \V nn + t nn , h s = h T = \{y nn - V nnn ). 
The first two terms in this effective pseudo-spin Hamil- 
tonian describe the intra-dimer interactions, while the 
last two ones describe its coupling with the off-domain- 
wall surroundings. It is worth noting that the condition 
hs = M means that the effective magnetic field produced 
by this surroundings turns into zero that provides the 
particle-hole symmetry of the dimer physics, in particu- 
lar, for the domain wall doping. Actually, we arrive at 
an effective singlet-triplet models Both S and T op- 
erators have a rather simple physical sense: the former 



corresponds to the total "quantum" charge of the dimer, 
while the latter does to the total "quantum" dipole mo- 
ment. Strictly speaking, the diagonal order parameters 
(S z ) and (T z ) describe the charge and dipole density, 
while the off-diagonal order parameters (S±) and (T±) 
describe the corresponding phase ordering. It is of pri- 
mary importance to note that these order parameters are 
not independent because the respective operators obey 
the kinematic constraint: 

S 2 + T 2 =3; (S-T) =0, 

stemmed from a simple spin algebra. It should be noted 
that there are two operators: T = si — S2 and J = [si X 82] 
that change the pseudo-spin multiplicity with their ma- 
trices being symmetric and antisymmetric, respectively: 

(00|T,|10) - (10|T 2 |00> = 1; 

(00|J,|10) = -(10|J,|00)=i. 

The Hamiltonian (|llfl points to the competition of S- and 
T-orders in the ring-shaped domain wall. 

FigEl shows a step-by-step formation of the energy 
spectrum of such a two-site cluster, or dimer. The ar- 
rows mark the dipole-allowed transitions that could be 
revealed in optical spectra. The tunnel states |00) and 
1 10) describing purely " half-filled" dimer states are mixed 
due to the electric field with the mixing level governed 
by the ratio \h,T/2t nn \. If | /it | 3> 2|i„„|, we arrive at 
classical Neel-like "up-and-down" dimer state. 

Hereafter, to describe our singlet-triplet quantum 
pseudo-spin system we start with trial functions 



ip = c 0oo + ^2( a j + ibj)ipj 



(12) 



where the spin functions \1M) in Cartesian basis are 



used: ij) z = 1 10 > and t/j. 



x.y 



|11 > ±|1 - 1 >)/y/2. 
The conventional spin operator is represented on this ba- 
sis by a simple matrix: 



(^ t \Sj\-ijjk) 



and for the order parameters one easily obtains: 

(S) = -2[a x b], ({SiSj}) = 2(% - a iaj - b^) (13) 

given the normalization constraint |c | 2 + a 2 + b 2 = 1. 
Thus, for the case of spin-1 system the order parameters 
are determined by two classical vectors (two real compo- 
nents of one complex vector c = a + ib from l|12[)). The 
two vectors are coupled, so the minimal number of dy- 
namic variables describing the S — 1 spin system appears 
to be equal to four. For the dimer as the singlet-triplet 
center we have additional unconventional T,J-order pa- 
rameters. For the respective averages we can easily ob- 
tain 



c c 



c c*). 



Given the real value of cq parameter we arrive at a simple 
form:(T) = 2c a, and (J) = -2c b. 
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FIG. 6: Pseudo-spin orientation in different symmetry sup 
components are shown. 

C. Pseudo-spin formalism and MFA description of 
the domain wall 

Now we may proceed with the inter-dimer coupling in 
the domain wall. Introducing the symmetrized superpo- 
sitions of S and T vectors, we can write the pseudo-spin 
Hamiltonian of the inter-dimer interaction as follows 

H dd = ^t nnn [S 2 (a lg )-S 2 (b lg ) + (S(e u )T(e u )) + T 2 (b 2g ) 

-T 2 (4,)] + \{Vnnn ~ t nnn )[S 2 (a lg ) - S*(b lg ) + 

(S z (e u )T z (e u )) + T 2 (b 2g ) - T 2 (b' 2g )], (14) 

where 

S(a lg ) = (Si + Sn + S ln + S IV ); 

S(6ig) = (Si — Sn + Sin — Siv); 

S(e u x) = (Si + Sn - Sm - S/y); 
S(e u y) = (-Si + Sn + Sm - S/y); 
T(b 2g ) = (Ti + T n - Tin - Tiv); 




for S- and T-order. Only the in-plane pseudospin 



T(&2 9 ) = (Tj - Tn - Tin + Tiv); 
T(e u x) = (Ti + Tn + T ln + T IV ); 

T(e u y) = (Ti - Tn + Tn, - T /y ); 

(see Fig0J). The chemical potential fj, for the system of 
four dimers in a domain wall is determined now by the 
condition: J2t=i(^i) ~ Ojil f° r undoped and singly- 
doped (boson/hole) domain wall, respectively. 

The respective mean values we can address to be or- 
der parameters that describe the subtle structure of a 
domain wall: as regards the diagonal order (S z (a\ g )) 
specifies the full charge; (S z (e u )) and (T z (e u )) do the 
electric dipole moments, (S z (b\ g )) the component of the 
quadrupole momentum, and one might introduce the 
higher order multipole moments. As regards the off- 
diagonal order we should, in general, proceed with the 
three types (ai g , e M , b\ g , or s,p,d) of the S-order, and 
three types (e u ,b2 g (b 2g ), or p, d) of the T-order, which 
can be defined as follows: 

(i_(7))=^e^,(f_(7))=p^. (15) 

Noteworthy to mention the kinematic constraint that 
couples different order parameters. In the absense of an 
external magnetic field the energy does not depend on b 
vector, it is restricted only to lie in xy-plane. In other 
words, we deal with the uncertainty of J-order parameter. 

The mean-field domain-wall ground state corresponds 
to the coexistent S(ai g ), T(t»2 ff ) modes given the negative 
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sign of the transfer integral t nnn or S(bi g ), T(& 2g ) modes 
given the positive sign of the transfer integral t nnn . The 
relation between S- and T-mode weight is specified by 
the relationship between the nn transfer integral t nn and 
electric field Iit- If we assume hx — 0, then given t nn > 
the dimers have purely S — 1 ground state, and we arrive 
at S-type off-diagonal order. In contrast, given t nn < 
we arrive at T-type off-diagonal order. In general, ne- 
glecting the ST-mixing term leads to two independent 
phase order parameters. Generally speaking, the inter- 
ference dipole-dipole ST-term would result in ST-mixing 
accompanied by the constraint on the phase order param- 
eters with an appearance of a noncollincarity effect. Thus 
we arrive at the conclusion that the symmetry of the or- 
der parameter distribution in the domain wall would be 
specified only by the sign of the transfer integral. In addi- 
tion, we see that the problem of the order parameter asso- 
ciated with our bubble domain is much more complicated 
than in conventional BCS-like approach^ due to its mul- 
ticomponent nature. Moreover, we deal with a system 
with different symmetry of low-lying excited states and 
competing order parameters that implies their possible 
ambiguous manifestation in either properties. The low- 
symmetry crystalline electric fields or crystal distortions 
would result in a mixing of the order parameters with 
different symmetry. The bubble domain (Fig0J yields a 
simple and instructive toy model to describe such effects. 

It is worth noting that the flux quantization effects 
for the bubble domain, in particular, the localizing ef- 
fect of the magnetic field on the moving bosons, are ex- 
pected to be observed only for rather large fields. Indeed, 
nanoscopic atomic systems such as a square plaquette 
with a size of around 10 A require a field of H 10 2 
Tesla for a half flux quantum per plaquette. 

Above we have concerned a simple one-center topo- 
logical defect. However, in general one has to make use 
of more complicated topological excitation like a multi- 
center BP skyrmion^i The question arises, whether such 
an entity may be described as a system of weakly cou- 
pled individual one-center defects? In practice, it seems 
to be a rather reasonable approach only for slight de- 
viation from half-filling, when the mean separation be- 
tween doped particles is much larger than the effective 
domain size. The bubble domain may be addressed to 
be well isolated only if its environment does not contain 
another domain(s) which could be involved to a "dan- 
gerous" nearest-neighbor inter-dimer coupling. Each do- 
main in FigQ]has z = 12 of such neighbors. Hence the 
concentration of well isolated domains can be written as 
follows: Po(x) = x(l — x) z , where z is a number of "dan- 
gerous" neighbors, x = An;, the boson concentration. 
The Pq (x) maximum is reached at xq = ^-j-j- . In our case 
xo = 1/13, or f» 0.077. With increasing doping the de- 
viation of Po(x) from the linear law rises. On the other 
hand, knowing the effective domain area Sd ~ 9a 2 we can 
roughly estimate the limiting concentration of the single- 
domain model description to be x max ~ 1/9, f=a 0.11. The 
inter-domain coupling includes both long-range SS- and 



short-range TT-terms. We should emphasize a strong 
anisotropy of this coupling. For nearest neighbors along 
[1,0] and [0,1] directions namely the short-range (dipole- 
dipole) TT-coupling would result in an effective suppres- 
sion of the main S ordering (see FigQJ , in contrast with 
[11] direction. 

Thus we may conclude that the above mentioned 
simple continuous model for Josephson junction arrays 
should be strongly modified to describe the multi-domain 
configurations in 2D quantum hc-BH lattices, both as re- 
gards the subtle internal domain structure, the competi- 
tion of different order parameter, and the anisotropy of 
Josephson coupling. 



V. CONCLUSIONS 

In conclusion, the boson/hole doping of the hard-core 
boson system away from half-filling is assumed to be a 
driving force for a nucleation of a multi-center skyrmion- 
like self-organized collective mode that resembles a sys- 
tem of CO bubble domains with a Bose superfluid and 
extra bosons both confined in domain walls. Such a topo- 
logical CO+BS phase separation, rather than an uniform 
mixed CO+BS supersolid phase, is believed to describe 
the evolution of hc-BH model away from half-filling. 
Starting from the classical model we predict the prop- 
erties of the respective quantum system. In frames of 
our scenario we may anticipate for the hc-BH model the 
emergence of an inhomogencous BS condensate for super- 
high temperatures T^ps < t, and 3D superconductivity 
for rather high temperatures T c < J < t. The system 
is believed to reveal many properties typical for granular 
superconductors, CDW materials, Wigner crystals, and 
multi-skyrmion system akin in a quantum Hall ferromag- 
netic state of a 2D electron gas. Topological inhomogene- 
ity is believed to be a generic property of 2D hard-core 
boson systems away from half-filling. Such a behavior 
represents a boson counterpart of the so-called topolog- 
ical doping being a general feature of Mott-insulator or 
2D fermion Hubbard modeliii 

Despite all shortcomings, MFA and continuous approx- 
imation are expected to provide a physically clear semi- 
quantitative picture of rather complex transformations 
taking place in bare CO system with doping, and can be 
instructive as a starting point to analyze possible scenar- 
ios. First of all, the MFA analysis allow us to consider 
the antiphase domain wall in CO phase to be a very ef- 
ficient ring-shaped potential well for the localization of 
a single extra boson (hole) thus forming a novel type of 
a topological defect with a single-charged domain wall. 
Such a defect can be addressed as a charged skyrmion- 
like quasiparticlc which energy can be approximated by 
its classical value for CO bubble domain. It is of great 
importance to note that domain wall simultaneously rep- 
resents a ring-shaped reservoir for Bose superfluid. 

Unfortunately, we have no experience to deal with 
multi-center skyrmions as regards its structure, energet- 



12 



ics, and stability. It should be noted that such a texture 
with strongly polarizable centers is believed to provide an 
effective screening of long-range boson-boson repulsion 
thus resulting in an additional self-stabilization. Nucle- 
ation of topological phase is likely to proceed in the way 
typical for the first order phase transitions. 

The role played by quantum effects and lattice dis- 
creteness has been illustrated in frames of the sim- 
plest nanoscopic counterpart of the bubble domain in 
a checkerboard CO phase of 2D hc-BH square lattice. 
It is shown that the relative magnitude and symmetry of 
multi-component order parameter are mainly determined 
by the sign of the nn and nnn transfer integrals. The 
topologically inhomogeneous phase of the hc-BH system 
away from the half-filling can exhibit the signatures both 
of s, d, and p symmetry of the off-diagonal order. The 
model allows us to study the subtle microscopic details of 



the order parameter distribution including its symmetry 
in a real rather than momentum space, though the prob- 
lem of the structure and stability of nanoscale domain 
configurations remains to be solved. The present pa- 
per establishes only the framework for analyzing the sub- 
tleties of the phase separation in a lattice hc-BH model 
away from half-filling. Much work remains to be done 
both in a macroscopic and microscopic approaches. 
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